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D . Abstract 



We combine the resummations for radiative corrections and for the heavy 
quark expansion to study the inclusive radiative decay B X^'y. The 
infrared renormalon ambiguity is also taken into account. Including both 
theoretical and experimental uncertainties, we determine the allowed domain 
for the HQET parameters A and Ai centered at A = 0.65 GeV and Ai = —0.71 
GeV2. 
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With the progress of the heavy quark effective theory (HQET), we have 
gained better insight to the dynamics of heavy systems. For example, the 
inclusive semileptonic B-meson decay rates can be expanded in powers of 
l/mb- To 0(l/m^) (ze within an accuracy of about 5%), three parameters 
A = mB—Trih, Ai = {{iDY) and A2 = {<y-G) are relevant. Except for A2 which 
can be obtained directly from the B*-B mass splitting, A and Ai are poorly 
determined. Recently, we have formulated the perturbative QCD (PQCD) 
approach to the inclusive semileptonic B-meson decays, which combines the 
resummation technique and the HQET based operator product expansion [Q]. 
In this approach the charged lepton spectrum is expressed as the convolution 
of a hard subprocess with a jet function and a universal i?-meson distribution 
function. 

In this letter we shall determine the parameters A and Ai, which are 
equivalent to the first and second moments of the 5-meson distribution func- 
tion, respectively, from the photon energy specrum of the radiative decay 
B —>■ Xg'-f p[. It should be emphasized that the relation A = — rrib, 
and thus the extraction of the pole mass nib of the b quark, suffer the in- 
frared (IR) renormalon ambiguity of power Aqcd/"^;) 0- Hence, we re- 
gard A as being related to the first moment of the distribution function, 
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instead of to the b quark mass. This viewpoint is satisfactory enough in 
the sense that the moments extracted here can be consistently employed to 
make predictions for other processes because of the universality of distri- 
bution functions. This makes possible a model-independent determination 
of the Cabibbo-Koboyashi-Maskawa (CKM) matrix elements \Vcb\ and \Vub\ 
from the inclusive semileptonic decays. 

On the other hand, the ambiguity of 0(Aqcd/'^6) in the definition of 
the pole mass turns out to be cancelled by the ambiguity contained in loop 
corrections , when one evaluates the total decay widths of heavy mesons. 
In the factorization theorem formulated at the meson level, we adopt the B- 
meson mass rriB and thus avoid the ambiguity in nib. It is then found that the 
IR renormalon ambiguity appearing in the perturbative resummation starts 
at 0{AQQj^/ml), consistent with the conclusion in [Q. This 0{AQfjY)/ml) 
ambiguity will be taken into account as a source of theoretical uncertainties. 
After including the uncertainty from the experimental data, we determine 
A = O.Q5t°ofo GeV and Ai = -0.7ll?:I[! GeV^. 

The first two moments of the photon energy spectrum of the decay 
B — > Xg'j have been computed using the expansions in l/rUb and in 
[^]. Without considering the errors of data, the parameter A ^ 450 MeV 
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was extracted. In our approach the PQCD expansion in ag is resummed 
up to next-to-leading logarithms, and the nonperturbative heavy quark ex- 
pansion in l/rrih is resummed into the 5-meson distribution function. Fur- 
thermore, both the IR renormalon ambiguity as a source of theoretical un- 
certainties and the errors of data as a source of experimental uncertainties 
are included. Therefore, our analysis is more complete. In |^ the contri- 
butions to the B Xg'-f decay from the b ^ sg transition through the 
s — > 7 and (? — > 7 fragmentation functions were studied, and found to be 
negligible. This observation hints that we concentrate only on the b s'j 
transition. Single logarithms in the fragmentation functions were summed 
using the renormalization-group (RG) method. In this work we employ the 
more sophiscated resummation technique to organize the double logarithms 
involved in the b ^ s'j decay. 

The one-loop constraint on A and Ai, A > [0.32-0.07(Ai/0.1GeV2)] GeV, 
has been obtained from the first moment of the invariant mass spectrum of 
the B —>■ XJu decay using the same expansion in both l/m^ and 0. This 
constraint, however, defines only an open domain. QCD sum rules are an 
alternative approach to the determination of the nonperturbative parameters, 
which give Ai = —0.6 ±0.1 GeV^ [0. It is also worthwhile to mention the 
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lattice extractions of rhblmb) = 4.17 ± 0.06 GeV and mb{mb) = 4.0 ±0.1 
GeV 1^, which are close to the lower bound of our predictions. It is obvious 
that our predictions are not only concrete, but consistent with the conclusions 
in the literature. 

The effective Hamiltonian for the process 6 — 57 written in terms of 
dimension-6 operators is 

H^sib - .7) = -^Xtj:C,{f,)0,{f,) , (1) 

with Gp the Fermi coupling constant, Cj(/i) the Wilson coefficients evaluated 
at the scale /i, and Aj = VtbV^* the product of the CKM matrix elements. 



The definition of the operators Oj is referred to |TI]]. Here we show only the 
relevant operator 



lOTT^ 



with R = {1 + 75)72 and L = (1 — 75)/2. Since the mixing between the 
operators Oi, ...,Oe and the operators Of and vanishes at one- loop level 
under the infinite renormalization, two- loop calculations of 0{eg'j.) and 0{gl) 
are necessary in order to obtain Cy{h) and Cs{fi) in the leading logarithmic 
approximation, where e and Qs are the electromagnetic and strong couplings, 



respectively. Recently, the 0{as) corrections to C7 and Cs have been com- 
puted in [O. For the b s'-f transition, the contributions from all the 
operatorrs can be included by simply employing the effective Wilson coeffi- 
cient C^^ = Cj + QdC^ + 3(5dC*6 EOl ) Qd being the charge of the d quark. 



The analysis in [|IOl involves only lowest-order PQCD corrections to the 
spectator model, and thus addresses nothing on the nonperturbative param- 
eters. However, the explicit 0{as) expressions for the decay width of 6 — 57 
help to understand our formalism. The infrared single pole appearing in the 
one-loop calculation is absorbed into the distribution function, the double 
logarithms In^ (ms/mt) with the s quark mass, are absorbed into the jet 
function, and the single logarithms ln(/i/mb) are absorbed into the hard part 
[0]. Employing the resummation technique and the RG method, these large 
corrections are grouped into a Sudakov factor. This systemmatic summation 
of large logarithms has not been achieved in the literature. The distribution 
function can be regarded as the consequence of the resummation of the heavy 



quark expansion |12 



The factorization formula is then written as 
1 dV 



r(o) dE^ 



niB J dz Jbdbf{z,b)J{z,x,mB,b)H{x)exp[—S{mB,b)], (3) 



with the tree-level decay rate T^^\ the jet function J, the hard part H, and 
the Sudakov exponent S given by 

r(°) = ^|Cf (m^)G^A,p«em , (4) 



J = Jo(V z — xniBh) , H = , (5) 
dp fP dfi 5 f^B dfias{fi) 



s=2r'j^r '-^AM,)) - r , (6) 

Jl/b P Jl/b 11 3 Jl/b IX 71 

(7) 



A = Cf— + 

IT 



/b p Jl/b /i 6 Jl/b fx TT 

67 TT^ 10 



9 3 27 

with Cp = 4/3 the color factor. The variable x is defined by x = 2E^/mB, 
being the photon energy. Note that the second term of A is a two-loop result 



[13]. For consistency, the two- loop expression of asifi) is inserted into the 
integral of S. Note that the 5-meson mass m^, instead of the 6-quark mass 
rrib, appears in the expression of T^^^ as stated before. The choice of the scale 
niB for Cj^, the same as the upper bound of the evolution of the Sudakov 
factor, follows the three-scale factorization theorem developed recently |]14 |. 



The argument z of the function f{z,b) is the momentum fraction, and b 
is the conjugate variable of the transverse momentum carried by the b quark. 
Including these transverse degrees of freedom, we set nis to zero, and let 1/6 
play the role of an infrared cutoff in the resummation for the jet function 
as shown in Eq. (|^). Since the intrinsic b dependence (the perturbative b 
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dependence has been collected into the exponent S) is not known yet, we 
take the ansatz f{z,b) = f{z) exp(— which leads to S ^ as 6 — 
according to the definition f{z,b = 0) = f{z). It is also natrual to assume 
S > for all b from the viewpoint that the b quark is bounded inside the 
B meson. Hence, the intrinsic b depedence provides further suppression. 
The nonperturbative function f{z), identified as the B-meson distribution 
function, can be expressed as the matrix element of the b quark fields, whose 
first three moments are |]T2| 

f{z)dz = 1 , 

Jo 

f f{z){l - z)dz =— + 0(A^cd/O , 
Jo niB 

f f{z){l - zfdz = -V - t) + ^(Aqcd/^I) • (8) 

JO ''''B \ / 

Though the exponent E is unknown, we can, however, extract its leading 
behavior by means of the IR renormalon analysis. Note that the pertur- 
bative Sudakov factor e""^ in Eq. (||) becomes unreliable as the transverse 
distance b approaches 1/Aqcd- Near this end point, as{l/b) diverges, and 
IR renormalon contributions are significant. We reexpress the RG result of 
the evolution of the distribution function, which is contained in the second 
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term of S, as 

W = exp 



47rCf 



(9) 



(27r)4 (t; . /)2 

The loop integral corresponds to the correction from a real soft gluon at- 
taching the two valence b quarks, whose propagators have been replaced 
by the eikonal lines in the direction v = (1,1,0) |[T^. The tensor N^" = 
g/iu ^j^/j-i^ _|_ [fj'fi'^'j . n'^l^^l^ / (n ■ /)^ comes from the gluon propagator 
in axial gauge n ■ A = Q. We have set the argument of the running to t^, 
which is conjugate to the scale h of the distribution function. 

Substituting the identity as(/|') = tt /q°° (icrexp[— 2cr/9i ln(Zr/AQCD)], A = 
(33— 2ri/)/12, into Eq. (P), and performing the loop integral forn oc (— 1, 1, 0) 



13 , we obtain 



W = exp 







2(T/3; 



(10) 



2 ; r(i + (T/3i; 

It is easy to observe that the pole of T{—aPi) at a ^ gives the perturbative 
anomalous dimension of the distribution function appearing in Eq. (^. The 
extra poles at a l//?i, 2//?i,..., then correspond to the IR renormalons, 
giving corrections of powers 6^, 6^,..., respectively. These renormalons gener- 
ate singularities, which must be compensated by the nonperturbative power 



correctons in order to have a well-defined perturbative expansion. Note that 



the renormalon ambiguity starts with the power (6Aqcd)^, instead of &Aqcd 
[0], because the function W does not include the self-energy corrections to 
the b quark, which vanish under the eikonal approximation. 



Using the "minimal" ansatz |jT5|], ie picking up only the leading (left- 
most) renormalon contribution, we parametrize the exponent by = 
c'm^fo^, corresponding to the fact that the power corrections start at 0{b^). 
Certainly, other parametrizations consistent with this fact are equally good. 
The IR renormalon ambiguity in other terms of S can be absorbed into 
Combined with the Sudakov exponent, which is approximated by S* ~ 
0.025m^6^, Eq. (J^) gives the differential branching ratio of the decay B — >■ 

dBR /-i 



dE'y C J X 



r— I rf2/(z)e-(^-")/^ , (11) 

with c = 4c' + 0.1 and r = 2r(°) rs/mB. By varying the parameter c under 
the constraint c > 0.1, the theoretical uncertainty arising from our formal- 
ism is taken into account. Assuming the values itlb = 5.279 GeV, Gp = 
1.16639 X 10-^ GeV"^ [|T§, Cf (rris) = -0.306 ± 0.050, |At| = 0.040 ±0.004, 



= 1/130 ||T0|, and the i?-meson lifetime tb = 1.60 ± 0.03 ps ||T7 



we 



have r = 1.90 x 10 GeV ^ with about 50% uncertainty. We find that the 
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allowed values of c run from 0.1 to 0.2, ie d associated with the nonpertur- 
bative exponent is a small number. It indicates that perturbative corrections 
are more important than nonperturbative corrections, consistent with the 
requirement from heavy quark symmetry. 

Since two parameters are to be determined, the function ji^z) is modeled 
by the two-parameter form proposed in |l|, 

/w - ^ [(---Iplp • '''' 

where is the normalization constant, and the parameters a and e will be 
adjusted to fit the CLEO data within errors. Note that the specific form 
of f{z) is not essential. Any two-parameter wave function serves the same 
purpose. A and Ai are then obtained from the moments of f{z) using Eq. (H). 
The numerical analysis proceeds in the following way: We choose the central 
value r = 1.90 x 10^^ GeV^^ first, and start with a value of a < 1. For an 
arbitrary value of e, we vary c under the constraint c > 0.1. If there exists 



a finite range of c such that the predictions from Eq. (11) fall into the error 
bars of the CLEO data 0, this e is acceptable. Repeating this procedure 
for different e, an allowed range of e is found. By changing a, we obtain the 
allowed domain for a and e, from which the corresponding domain of A and 
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Ai is determined. 

Results are shown in Fig. 1, in which the sohd curves correspond to 
a = 0.91, 0.92,..., and 1.00 from right to left, and a point on each curve 
corresponds to a value of e. For a < 0.90, no acceptable e exists. These 
curves form the corresponding allowed domain of A and Ai. Including the 
errors of r, the domain enlarges by 20%. If there is no the constraint c > 0.1, 
the upper bound of A will increase 30%. The dashed curves are quoted from 
[^], with the left one and the right one obtained from the first and the second 
moments of the invariant mass spectrum of the B —>■ XJiy decay, respectively. 
The space above the dashed curves is the allowed domain in 0. The third 
information from the ratio of partial widths Rr = T{B XcTv) /T{B — >■ 
X^ev) introduces a constraint from the top. However, this constraint is not 
yet convincing [§, and still leaves the allowed domain an open one. Though 
the overlap of our results with those in is not very large, it is reasonable 
to claim that they are consistent with each other, because the theoretical 
uncertainties of the approach in was not estimated. 

We take the middle of the a = 0.95 curve as the central values, which 
lead to A = 0.65to,tl GeV and Ai = -O.ntiH GeV^. The bounds of these 
extractions are indeed very large. However, we emphasize that they are the 
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consequence of the inclusion of as more as possible theoretical and exper- 
imental uncertainties into our formalism. These bounds will definitely be 
reduced when, for example, more accurate data are available. If the errors of 
the data become half of current ones, we shall obtain A = 0.65^0^3 GeV and 
Ai = -0.71 

-im GeV^. As a simple estimation, we extract from A the 6-quark 
mass nib = 4.63;j!o;3o GeV, which certainly suffers the IR renormalon ambi- 
guity. Employing the central values, we determine the i?-meson distribution 
function 



which can be used to determine the inclusive semileptonic and nonleptonic 
decay spectra in the future. 

To highlight the resummation effect, we present the predictions for the 
photon energy spectra of the B Xg'y decay derived from the nonperturba- 
tive HQET distribution function alone, and from our formula including the 
perturbative resummation in Eq. (|TT]). The expression of the former is sim- 
ply dBR/dE^ = rx^/(x). Substituting r = 1.90 and Eq. ([TBp into the above 
expression and into Eq. ([TTl) with c = 0.15, we obtain the spectra shown in 
Fig. 2. The CLEO data are also displayed. It is obvious that the spectrum 




(13) 



[(2- 0.95)2 + 0.0034z] 



2 ' 
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from f{x) has a sharp peak near the high end of Ej, which satisfies the 
HQET, but is in conflict with the data. The predictions match the data only 
after the suppression effect is included, which possess a softer profile. Note 
that a naive fitting without the suppression effect needs a board distribution 
function, leading to a value of nib as small as 3.0 GeV. 

The branching ratio can be evaluated simply by integrating Eq. ([iTD over 
E^. We obtain BR = 2.80T 

oio ^ 10 where the central value corresponds 
to the solid curve in Fig. 2 (c = 0.15), and the upper bound and lower 
bound to c = 0.13 and c = 0.25, respectively. Our prediction is close to the 
standard-model estimation (2.8 ± 0.8) x 10~^ [18|, and to the CLEO data 
(2.75 ± 0.67) X 10~^ from the 5-reconstruction analysis 0. 

This work was supported by the National Science Council of the Republic 
of China under Grand Nos. NSC-85-2112-M-194-009 and NSC-85-2112-001- 
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Figure Captions 



Fig.l: The allowed domain for A and Ai. The dashed curves are the con- 
straints quoted from 0. 

Fig. 2: The photon energy spectra derived from the 5-meson distribution 
function alone (dashed curve), and from the inclusion of the suppression 
effect (solid line). The CLEO data are also shown P|. 
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